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Conquer the fine structure splitting of excitons in self-assembled InAs/GaAs 
quantum dots via combined stresses 
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Eliminating the fine structure splitting (FSS) of excitons in self-assembled quantum dots (QDs) is essential to 
the generation of high quality entangled photon pairs. It has been shown that the FSS has a lower bound under 
uniaxial stress. In this letter, we show that the FSS of excitons in a general self-assembled InGaAs/GaAs QD 
can be fully suppressed via combined stresses along the [110] and [010] directions. The result is confirmed 
by atomic empirical pseudopotential calculations. For all the QDs we studied, the FSS can be tuned to be 
vanishingly small (< 0.1 fJ-eV), which is sufficient small for high quality entangled photon emission. 



There has been continuous interest in finding an effi- 
cient entangled photon source for quantum information 
applications. Benson et al. proposed that a biexciton 
cascade process in a self-assembled quantum dot (QD) 
can be used to generate the "event-ready" entangled pho- 
ton pairs. 1 In this scheme, a biexciton decays into two 
photons via two paths of different polarizations \H) and 
\V). If the two paths are indistinguishable, the two emit- 
ted photons are polarization entangled. 1 ' 2 Unfortunately, 
the | if)- and |V)-polarized photons may have a small 
energy difference, known as the fine structure splitting 
(FSS), 3,4 due to the asymmetric electron-hole exchange 
interaction in the QDs. 5 ' 6 The FSS is typically about -40 
~ +80 [ieV in the InAs/GaAs QDs, 7 which is much larger 
than the radiative lincwidth (~ 1.0 ficV), and therefore 
provides "which way" information about the photon de- 
cay path that destroys the photon entanglement. 2,8 Al- 
though, it is possible to "cherry-pick" a QD that has 
tiny FSS from a large amount of QDs, 2 it is highly de- 
sired that the FSS can be tuned by the external fields in 
a controlled way. Indeed, it has been successfully demon- 
strated that the FSS can be tuned by electric fields, 9-13 
magnetic fields, 14 and uniaxial stress. 15 ' 16 

Applying a stress is an effective way to reduce the FSS 
in a QD experimentally. 15,16 For an ideal QD with Gi v 
symmetry, the FSS can be reduced to exact zero when 
the stress is applied along the [110] direction. 17 However, 
for a general QD with C\ symmetry, there is a lower 
bound for the FSS 17 because the two bright states be- 
long to the same symmetry representation. Gong et al. 
derived a general relation between the FSS, the exciton 
polarization angle and the uniaxial stress. They have 
shown that the FSS lower bound can be predicted by the 
polarization angle and the FSS under zero stress. Similar 
relationship is also found when the FSS is tuned by verti- 
cal electric field, 13 where the minimal FSS under electric 
field can be as small as 0.7 fieV. However, for most of 
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FIG. 1. The InGaAs QD is embedded in a 60 x 60 x 60 GaAs 
matrix. External stresses are applied along the [010] and [110] 
directions simultaneously. 



the dots, the minimal FSS under vertical electric field is 
still larger than the required splitting to observe entan- 
gled photons emission (see Supplementary Information 
in Ref. 13, where 16 of the 22 studied dots have the min- 
imum FSS above 2.5 /xeV). One may ask whether we can 
further reduce the FSS of a general QD by the external 
electric fields or stresses, to nearly zero. 

In this Letter, we demonstrate that indeed the FSS of a 
general self-assembled InGaAs/GaAs QDs with C\ sym- 
metry can be tuned to be vanishingly small by applying 
combined stresses along the [110] and [010] directions. 
The results are confirmed by atomistic empirical pseu- 
dopotential calculations. 19,20 The minimal FSS is gener- 
ally below 0.1 ^eV, which is sufficient for high quality 
on-demand entangled photon emissions. 

In a previous work, Gong et al. proposed an exciton 
Hamiltonian under uniaxial stress, 18 



H(n,p)=H 2v + V 1 +V s (n)p, 



(1) 



a ' Electronic mail: hclx@ustc.edu.cn 



where n is the external stress direction, and p is the mag- 
nitude of the stress. H^v represent the Hamiltonian of an 
ideal QD with C^v symmetry, whereas V\ lower the dot 
symmetry to C±, due to local structure deformations, al- 
loy distribution 21 and interfacial effects 22 etc. V s (n)p is 
the potential change due to the external stress. Note 
that the lattice deformation under external stress (< 10 2 
MPa) is generally smaller then 0.1%, thus in the above 
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Hamiltonian the term 0(p 2 ) is omitted. The eigenvec- 
tors of the two bright states of H 2v are |3) = |r 2 ) + 
and |4) = ^2) — i\^4), with corresponding eigenvalues 
E3 and £4, respectively. In the absence of an in-plane 
magnetic field, the coupling between the dark states and 
bright states is negligible. We therefore write the Hamil- 
tonian in the space spanned by the two bright states of 



H = 



E 



-f 5 + a 3 p 
K + (3p 



E 



k + f3p 
— 6 + ot\p 



(2) 



where E+S = (3|ff 2o +Vi|3), E-S = (4| J ff 2 „+^i|4). a t = 
(i\V„(n)\i) (i=3, 4), K = (3|Vi|4) and p = (3|V s (n)|4). 
Because the Hamiltonian has a time-reversal symmetry, 
all parameters can therefore be set to real values for sim- 
plicity. We define a = a 3 — and 7 = a 3 + 0:4. The 
cigenstates of the two bright states are the linear com- 
bination of states, |3), |4) i.e., \ijj-) = cos#|3) + sin#|4) 
and = -sin6>|3) + cos(9|4). The FSS A{p) for QDs 
under uniaxial stress p is, 



A(p) = y/Mj3p + k) 2 + {ap + 28f 



(3) 



According to the symmetry analysis (and verified by the 
atomistic theory), one has (3 = when the stress is ap- 
plied along the [110] ([110]) direction, and a=0, if the 
stress is applied along the [010] ([100]) direction. 18 Since 
the uniaxial stress can not make (J3p + k) and (ap + 25) 
equal to zero simultaneously, there is always a lower 
bound for FSS when kS ^ 0, which is determined by 
2\k\ or 2\S\ depending on the direction of the stress. 

The above results suggest that diagonal and off- 
diagonal terms of the Hamiltonian matrix Eq. (2) can 
be tuned by the stress along the [110] and [010] direc- 
tions independently. Therefore, it is possible to further 
reduce the FSS by applying the stresses along the [110] 
and [010] directions simultaneously. To show this, we 
generalize the 2x2 model exciton Hamiltonian under uni- 
axial stress to that under the combined stresses, 



H = 



E + 5 + a 3 j)[ii ] k + Pp[ow] 

« + #P[010] E-S + O4P[ii ] 



The FSS for QDs under the combined stress is 



A(p) = ^/4(/3p [010 ] + k) 2 + (ap[uo] + W 



(4) 



(5) 



The polarization angle 
culated from, 

tan(0±) = - 



vs 



P[110] 



and 



P[oio] 



can be cal- 



-25 - ap[i 10 ] T A(p) 
2(/3p[oio] + K ) 



(6) 



Naively, one may expect that by choosing proio] = — K>/fl 
and proiol = —28 /a, the FSS can be tuned to be exactly 
zero. However, it should be note that the exact zero 
FSS is forbidden by the symmetry restriction for a QD 
with C\ symmetry. In the above model, we neglect the 
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FIG. 2. Change of FSS (a) and polarization angle (b) under 
the combined stress along the [110] and [010] directions for 
the lens-shaped Ino.6Gao.4As/GaAs dot. The scatters are the 
results of numerical calculations whereas the solid lines are 
the results of the two-level model. 



higher order terms, which may become important when 
FSS becomes very small. To see what is the minimal FSS 
one can obtain under the combines stresses, we perform 
atomistic empirical pscudopotential calculations 19,20 on 
the FSS under the combined stressed along the [010] and 
[110] directions. Remarkably, we find that the FSS of a 
general InGaAs/GaAs QD can be reduced to be vanish- 
ingly small (<0.1 fieV) under the combined stresses. 

The InGaAs/GaAs QDs are modeled by embedding 
the InGaAs dots in a 60 x 60 x 60 8-atom GaAs super- 
cell. The QDs are assumed to be grown along the [001] 
direction, on the top of the one monolayer InAs wetting 
layers. We apply combined stresses along the [110] and 
the [010] directions as shown in Fig. 1. The stresses can 
be either compressive or tensile depending on the struc- 
ture of QDs. The resultant deformation of the supercell 
under stresses is approximated by the linear superposi- 
tion of the deformation of the cell under individual uniax- 
ial stresses. Optimal atomic positions under stresses are 
obatained with valance force field method. 23 ' 24 We solve 
the Schrodingcr equation to obtain the single-particle en- 
ergy levels and wavefunctions using a strained linear com- 
bination of Bloch band method. 25 The exciton energies 
are then calculated via the many-particle configuration 
interaction method, 26 in which the exciton wave func- 
tions are expanded in Slater determinants constructed 
from all confined electron and hole single-particle states. 

We have calculated 8 InAs/GaAs QDs with different 
geometries, including (elongated-)lens, pyramid, trun- 
cated cone, etc. and different alloy compositions. Here 
we take a lens-shaped In^Gai-^As/GaAs dot with x = 
0.6, diameter D = 25 nm, and height h= 3.5 as an exam- 
ple. We obtain similar results for all other QDs. Fig. 2 
depicts the behavior of the FSS and polarization angle of 
the dot under combined [110] and [010] stresses. We plot 
the FSS and exciton polarization angles as functions of 
stress along the [110] direction under P[oin]=0, 20, 40 and 
54 MPa. For a given stress along the [010] direction, the 
FSS goes through an anti-crossing as a function of stress 
along the [110] direction, which is very similar to the re- 
sults under uniaxial stress. 18 The critical stress along the 
[110] direction changes little with the stress applied along 
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TABLE I. Parameters for In^Gai-xAs/GaAs QDs with different shape and compositions under combined stresses along the 
[110] and [010] directions. The unit of dot diameter D and hight h is nm. The fine structure A(p) is in unit of fieV and the 
combined critical stresses p c = (p[no],P[oio]) are m un it of MPa. 







Pesudopotential Calculation 


Model Parameters and Predictions 


QDs 




A(p c ) 


Pc 


a 
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P 




Pc 


Cone (x = 0.6) 
D = 25, h = 3.5 




0.084 


(58.0, -36.0) 


-0.15 


4.36 


0.06 


2.04 


(58.6, -36.1) 


Truncated Cone (x = 


= 0.6) 
















D B = 25 




0.032 


(11.3, -45.7) 


-0.14 


0.79 


0.05 


2.23 


(11.5, -46.0) 


D T = 20, h = 3.5 


















Elongated Lens (x = 


0.6) 
















D im = 20 




0.045 


(75.8, -19.8) 


-0.17 


6.37 


0.06 


1.24 


(76.0, -20.1) 


D fl ioi = 25, h = 3.5 


















Elongated Lens (x = 


0.6) 
















D lm = 25 




0.076 


(-36.8, -10.7) 


-0.17 


3.03 


0.06 


0.70 


(-36.6, -11.4) 


•D[iioi = 20, A = 3.5 


















Lens (x = 0.8) 
£> = 25 h = 3.5 




0.048 


(35.6, 14.0) 


-0.19 


3.46 


0.06 


-0.89 


(35.8, 14.1) 


Lens [x = 0.6) 
D = 25 h = 3.5 




0.039 


(58.5, 54.0) 


-0.14 


4.11 


0.05 


-2.79 


(58.5, 54.3) 


Lens (a; = 0.6) 
D = 25 h = 5.0 




0.083 


(-12.4, -43.6) 


-0.13 


-0.79 


0.045 


1.99 


(-12.2, -43.8) 


Pyramid (a; = 0.6) 
D = 25 h = 3.5 




0.075 


(27.5, 1.6) 


-0.13 


1.77 


0.051 


-0.092 


(28.0, 1.78) 



the [010] direction. However the minimal FSS changes 
dramatically for different stresses applied along the [010] 
direction, as we expected. When no stress is applied, the 
FSS of the QD is approximately 9.8 /zeV. When we ap- 
ply the stress along the [110] direction, one can reduce 
the FSS to approximately 5.6 /xeV, and 8.2 /ieV if the 
stress is applied along the [010] direction. In cither case, 
the FSS is much larger than 1 ficV to suitable for high 
quality entangled photon emitters. However, when com- 
bined stresses along the [110] and [010] directions are ap- 
plied, the lowest FSS we obtain is 0.04 /xeV at p[ 110 ]=58.5 
MPa and prrjio] =54 MPa. The polarization angles rotate 
clockwise or anticlockwise depending on the magnitude 
of stress along the [010] direction and the polarization 
angle at the critical stress always tends to be —45° or 
45°. The calculated FSS and exciton polarization angles 
are in excellent agreement with model predictions. The 
minimal FSS of all other QDs are also <0.1 fieV under 
the combined stresses (See Table I for critical stresses and 
the FSS minimum of all dots we studied). The tunable 
area in which FSS is less than 1 /xeV is proportional to 
2/a and In typical InGaAs/GaAs dots, \a\ ~ 0.1-0.4 
/ieV/MPa and \/3\ - 0.04-0.1 ^cV/MPa, 18 therefore the 
tunable area for FSS less than 1 [ieV is approximately 
5 x 10 MPa 2 , which is easy to control in experiments. 

When the stresses along the two directions approach 
the critical point at the same time, the FSS becomes very 
small, the higher order terms in Eq. (4) may not be ig- 
nored, and therefore the FSS can never be exactly zero. 
Such higher order terms are fully included in the atom- 
istic calculation. The comparison between the model pre- 
dictions and exact numerical simulations can reveal the 
role of these terms. As shown in Fig. 2, the numerical 



results (scatters) and two-level model predictions (solid 
lines), using the parameters fitted from the results of uni- 
axial stress, are in excellent agreement. This suggests 
that the higher order terms are indeed unimportant ex- 
cept extremely close to the critical stresses. 




Stress Along [1 1 0] (MPa) Stress Along [110] (MPa) 



FIG. 3. Left panel: The change of the total dipole moment 
under the stresses for the two bright states: \4>-) ( a ) an d \ip+) 

(c) , in the lens-shaped Ino.6Gao.4As/GaAs dot. Right panel: 
The change of the z-component of dipole moment under com- 
bined stresses for the two bright states: \tp~) (b) and |^>+) 

(d) . 

Figure 3(a), 3(c) depict the total dipole moment 
\p\ 2 =Px+pl+pl of the two bright states under com- 
bined stresses. The total dipole moments of the two 
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states change a little in opposite directions under the 
stresses, and the lifetimes of the excitons change accord- 
ingly. However, the sum of the total dipolc moments of 
the two bright states shows very little change (< 1%) for 
stresses up to 100 MPa. Interestingly, for applied stress 
P[oio] = 54 MPa (critical stress along the [010] direc- 
tion), the total dipole moments of the two states change 
suddenly near the critical stress pmo] = 58 Mpa. This 
is because, at P[oio] = 54 MPa, the dipole moments of 
the two states are approximately aligned along the [110] 
and [110] directions, with one increases whereas one de- 
creases under the stress along the [110] direction. When 
the stress approach the critical stress along the [110] di- 
rection, the polarization angle rotates very fast from 0° 
to —90°, i.e., the two states switch polarization, which 
causes the jump of the total dipole moments of the two 
states. There are also very small z components of the 
dipole moment p z because of the valence band mixing, 
as shown in Fig. 3(b), 3(d). Therefore, when we project 
the dipole moment into the x-y plane, the polarization of 
the two bright states are not exactly orthogonal to each 
other. The deviation is approximately 1° - 2°. 

Bennett et al. have successfully tuned the FSS by ap- 
plying a vertical electric field. 13 According to the sym- 
metry analysis, the effect of a vertical electric field is 
equivalent to the uniaxial stress along the [110] direc- 
tion. Therefore, it is also possible to tunc the FSS to 
approximately zero by applying a combined vertical elec- 
tric field and an uniaxial stress along the [010] direction, 
which might be more feasible experimentally We leave 
this for future study. 

To summarize, we have demonstrated that the FSS in 
general self-assembled QDs can be reduced to nearly zero 
under combined stresses along the [110] and [010] direc- 
tions. A wide range of stresses can reduce the FSS to less 
than 1 /xeV and thus provide great merit in experiments. 
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